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Abstract 

m 

■ The Dyson Brownian motion model for transitions to the CUE is considered. For initial eigenvalue 

probability density functions corresponding to the COE and CSE, the density-density correlation 

G\ ■ 

function between an eigenvalue at position A' initially (r = 0) and an eigenvalue at position A for 
general r, is calculated. Theoretical predictions for the asymptotic behaviour are verified, and the 
distribution for the initial rate of change of the eigenvalues is calculated. Also, for initial conditions 
corresponding to equi-spaced eigenvalues, the n-point equal parameter distribution is computed. The 
extension of the calculation of these quantities to a multiparameter version of the theory is given. 
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1 Introduction 

A problem of fundamental importance in studying quantum aspects of classically chaotic systems is the 
statistical properties of the energy spectra. Since the pioneering work of Wigner and Dyson (see e.g. [1]) 
it has been realized that the statistical properties are determined by the presence, or not, of a global time 
reversal symmetry T. Details of the matrix elements of the Hamiltonian are unimportant, apart from 
the global constraint imposed by the symmetry T . Thus the discrete portion of the energy spectrum 
can be studied by modelling the Hamiltonian by a random Hermitian matrix if there is no time reversal 
symmetry, and in the case of a time reversal symmetry by a random real symmetric matrix if T 2 = 1 
and a random real quaternion matrix if T 2 = —1. 

More complicated situations can arise. In particular the quantum system may contain a parameter 
which can continuously vary the spectrum within a universality class or indeed between universality 
classes (see e.g. [2,3]). For the latter situation, by considering parameter (r) dependent Gaussian random 
matrices (see e.g. [4]), it is possible to show that the eigenvalue probability density function (p.d.f.) 
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p(Ai, . . . , Ajv; t) satisfies the Fokker-Planck equation 



dp „ , d fdW n , d \ . ,, 



where 

N 



W = - J2 loglA.-Afel + i^Al. (1.2) 



2 

l<j<k<N j=l 

Here 7 provides a scale for r/Xj and /3 = 1, 2 or 4 depending on whether the transition is to orthogonal, 
unitary or symplcctic symmetry. With r = 0, p is specified as the eigenvalue p.d.f. of one of the three 
Gaussian random matrix ensembles, with a symmetry different to that after the transition. In the random 
matrix context, the Fokker-Planck equation (1.1) is referred to as the Dyson Brownian motion model. 

The objective of this paper is to provide the exact evaluation of some two-point correlation functions 
for the transition to unitary symmetry In Section 2 we revise exact results on this problem which are 
already known [5,6], as well as a theoretical prediction for asymptotic behaviours, the implication this 
prediction has for the integrated correlator of the variance of a linear statistic [7], and the calculation 
of the distribution function for the initial rate of change of the eigenvalues. In Section 3 we provide the 
exact evaluation for the density-density correlation between an eigenvalue at r = and an eigenvalue at 
general r, when the eigenvalue p.d.f at r = corresponds to orthogonal symmetry, while in Section 4 
the calculation will be repeated with the r = eigenvalue p.d.f. corresponding to symplectic symmetry. 
In both cases the expected asymptotic behaviour is verified, and the distribution function for the initial 
rate of change of the eigenvalues is calculated. In Section 5 we calculate the equal-parameter n-point 
correlations when the initial conditions correspond to equispaced eigenvalues, while in Section 6 we show 
how the results of the previous sections can be extended to apply to the multiple parameter case. 

2 Summary of known results 

2.1 Equal parameter distributions 

The so called equal-parameter distribution functions 

p n (X 1 ,...,\ N ;r):=N(N-l)...(N-n)( 11 f° dX^\ p(X u . . . , X N ; r) (2.1) 

Vz=™+i"'- 00 / 

have been calculated exactly for the Fokker-Planck equation (1.1) with (3 = 2, and initial conditions 

p(A 1 ,...,A JV ;0) = -^- e -^ 1 -^) (2.2a) 

^1N 

N/2 

p(X u ...,X N] 0) = - e -4W(x 2 M,...,x N ) TT 6{X2i _ X2ii) (Areven) (2.26) 

G 4(JV/2) fJl 

in refs. [5,6]. 

We are particularly interested in these distribution functions in the thermodynamic limit, when TV — > 
00 while the density of eigenvalues is fixed at some non-zero constant p. In the TV — > 00 limit, the density 
implied by (1.1) is (see e.g. [4]) the Wigner semi-circle 

2V2TV 



p(X) ~ — y/1 - A 2 /2TV, 
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so to obtain a bulk density p about the origin it is necessary to scale the eigenvalues Xj i— > irpXj/\ / 2N. 
This can be achieved by instead scaling 7 so that 



7 1 ► 2Nj/(np) 2 . 



(2.3) 



With this scale, the thermodynamic limit of the equal time distributions (2.1) exist and can be computed 
exactly for the initial conditions (2.2). In particular, for the initial conditions (2.2), the two-point 
distribution function is given by [5,6] (for consistency with our notation the parameter p 2 in [5] is identified 
with p 2 r/2j and r replaced by pX) 



P(2)(Ai, A 2 ) 



o(A) og(A,r) 



4l{\r) a(X) 



ag(A,r) 



a(A) — "i2 
-4{(\t) a(A) 



= ^-^(a(A) 2 -a*(A,T)a*(A,r)) 



(2.4) 



where 



A := Ai — A 2 a(A) := 



sin 7rpA 
npX 



and for the initial condition (2.2a) ([3 = 1) 

a 12( A ' T ) : = _7r / 

Jo 

a 2i(A,r) :: 
while for the initial condition (2.2b) 

«i 2 (A,r) := 



duusmirupXe* 2 " 2 ™ 2 /"' 
du 



shiTrupA 2 p 2 TU 2 h 
u 



«2l(A,T) 



(iwii sin 7rup\e~ 
1 f\ u ^up\^ p ^ h 

V J U 



(2.5a) 
(2.56) 

(2.6a) 
(2.66) 



2.2 Density-density and current-current correlations 

To define the density-density correlation, S(X, A';r) say, between an eigenvalue A' at r = and an 
eigenvalue A at general r, let 

G(A' 1; . . . , A' w |Ai, . . . , X n ;t) 

denote the Green's function for the Fokker-Planck equation (1.1) i.e. the solution which satisfies the initial 
condition 



JV 



G(X[, . . . , A'jv|Ai, . . .,X N :0) = JJ<5(A ; - Aj), 



(2.7) 



and let po(X[, . . . , A^) denote the initial eigenvalue p.d.f., and p^ (A, r) denote the one-point distribution 
(the density) defined by (2.1). Then 

N /.no N 



S(A,A';r) := (f[ / dAj) ^J(Aj - A')po(Ai, . . . , A^) 
i=i J -°° j=i 

" ,00 JV 

x II / rfA E «*( A J - A ) G ( A 'i > • • • ' Ajv|Ai, • • . , Ajv;t) - (A', 0)p (1) (A, r). 



(2.8 
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In the thermodynamic limit, with the scaling (2.3), 

P(i)( x , T ) = P 

for all r and we expect S to depend only on the difference A — A', so we write 

5(A,A';r) := S{X-X';t). 

On the basis of a non-local diffusion equation approximation to (1.1) Beenakker and Rejaei [7] have 
obtained a prediction for the small-A; behaviour of the Fourier transform 

/oc 
S(X- X';T)e lk(x - x 'U\. 
-oc 

This prediction states that 

S(k;r) ~ S(fc;0)e~" Vr|fe|/7 as k -» 0. (2.9) 
In particular, if the initial p.d.f. is 

' e-f> w (2.10) 



CfiN 

where W is given by (1.2), then [1] 

S*(fc;0) ~ ^4 as l fc l -» °> 
so taking the inverse transform of (2.9) gives 

S(X-X';r) ~ -l^^log((^r/ 7 ) 2 + (A-A') 2 ) for |A-A'|-oo. (2.11) 

When the initial eigenvalue p.d.f. is (2.10) with (3 = 2 (the Fokker-Planck equation then describes a 
perturbation of the unitary ensemble), the density-density correlation is well known (see e.g. [7]): 

S{X- X'-t) = p 2 J d Ul cosirmpXe 712 " 2 ™!/ 2 " 1 J du 2 cos iru 2 pX e -*V™2/27. (2.12) 

(with the eigenvalues A and A' chosen in the neighbourhood of the edge of the Wigner semi-circle, the 
exact evaluation of S is also known [8]). The asymptotic expansion of (2.12) agrees with (2.11). 

In this paper we will provide the exact evaluation of S(X — A'; r) for the initial conditions (2.2). For 
the initial conditions (2.2a) the above theory gives 

5(fc;r) ~ iM e --V|fc|r/7 as |Asj — ► 0, (2.13a) 
while for the initial conditon (2.2b) 

S(k;0) ~ 2^ as |fc| -» 0, 
(the factor of 2 accounts for the double degeneracy in (2.2b)), and so (2.9) gives 



J*i e -"V 

2lT 

Now consider the correlator R(t; a) of a linear statistic 



S(k;r) ~ U e -- 2 p 2 |fe|r/7 as |As| — ► 0. (2.136) 



N 

A(r;a) := ^ a(Xj(T)/a). 

3 = 1 
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The correlator is defined as 

R(t; a) := ((A(0; a) - (A{0; a)))(A{ T ; a) - (A(t; a))) 

/oo POO 
dX / d\'a(\/a)a(\'/a)S(\- A';r), (2.14) 
-oo J -co 

where the averages in the first line above are over both the eigenvalues at r = and the eigenvalues 
at general r. Beenakker and Rejaei [7] used the asymptotic formula (2.9) to compute the integrated 
correlator 

/•OO 

XA ■= / dxR(x 2 ;a), 
Jo 

for the initial p.d.f. (2.10), in the limit a — > oo (the quantity x 2 = r is identified in [7] as an external 
perturbation parameter). Thus, introducing Fourier transforms, we have 

POO POO 

Xa = tt dx dk\a(k)\ 2 S(k/a;x 2 ). 

271" Jo J -co 

In the limit a — > oo, use of (2.9) gives, after carrying out the integration over x, (eq. (4.22) of [7] with 

7 i ^ 7/Vo) 

1/2 ,.00 



1 



p OO 

/ dfc|a(/fc)| 2 £; 1/2 . 
Jo 



Another consequence of (2.9) is a formula for the current-current correlation, defined formally as 

C7(A,A';r) :={£ ^ ^ (A - A,(0))*(A - A fc (r))) 
j',fe=i 

In a Brownian motion model the derivatives cannot be computed directly. However, by considering the 
correlator (2.14) with A replaced by dA/dx, it can be shown that the current-current correlation is related 
to the density-density correlation by [7, eq. 4.14] 

C(k;x 2 ) = ~^S(k;x 2 ). (2.15) 

Thus, from (2.9), for small- 1 A; | and the initial condition (2.10), one obtains [7] 



C(k;x 2 ) ~ |^(l-2^ 2 p 2 |fc|x 2 /7)e- 7rV|fc|x2/7 , 
P7 



or equivalcntly 



C{\\':x 2 ):=C(\-\'-x 2 ) ~ log (W/ 7 ) 2 + (A - A') 2 ) (2.16) 



for large- 1 A — A'|. 

2.3 The distribution of the initial rate of change of the eigenvalues 

As commented in the above section, the parameter x :— r 1 ' 2 can be identified with an external perturbing 
parameter. For the initial conditions (2.2b) a quantity of interest [9] associated with this interpretation 
is the distribution of the scaled repulsion between degenerate levels as the perturbation is initialized: 



A(a) := lim lim ( — "S^ 5 



2 W a X 2j (T) - X 2j ^(r) 
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N/2 

(r))/r^) 



I 1 2 f 
= lim lim ( V / cW^-f^M-A^ 

r^OA^oo \27TiV ^7_oo 
I />oo /•oo 

= lim — / dpe*** / dApf 2) (A;r) e 
r^o 2np J_ OQ J_ OQ y ' 



lim — pj^ar^r) (2.17) 



P 

The exact result (2.4) gives [9] 

A(a) = (^) lf \ a2e ^ a2/2 - ( 2 - 18 ) 

A quantity of interest for general initial conditions is the distribution of the initial rate of change of 
the eigenvalues with respect to x: 

iXp/r 1 ' 2 



-1 />OG /'OC 

lim / dpe lpv dXS(X,T)e~ 



T-l/2 

= lim S^V) (2.19) 

T— >0 p 

Using the exact result (2.4) for initial conditions with unitary symmetry, (2.19) gives [10] the Gaussian 
distribution 

f( v )=ay /2e ~ v2j/2 - ( 2 - 2 °) 



3 Exact solution for initial conditions with orthogonal symme- 
try 

3.1 Green's function for /3 — 2 in periodic model 

Our interest is in calculating the correlation (2.8) in the thermodynamic limit. For this purpose we can 
replace the potential W in (1.2) by its periodic version 

W = - log|sin7r(A fe - Aj)/L|, -L/2 < Xj < L/2. (3.1) 

l<j<k<N 

This gives the Fokker-Planck equation of the Dyson Brownian motion model for the circular ensemble of 
unitary random matrices [11]. 

In general the equilibrium solution of (1.1) is 

P= 7 J-e-? w , (3.2) 

and the operator 

e PW/2£ e ~pW/2 

is Hermitian [12]. With (3 = 2 and W given by (3.1), a standard calculation gives 

1 N Fft 



G 



where 

Kn = I 

12 

The well known reason for the tractibility of the case (3 = 2 is thus revealed: the Hermitian operator on the 
r.h.s. of (3.3) is the Schrodinger operator for non-interacting particles. To reclaim the correct equilibrium 
solution (3.2) the particles are assumed to be impenetrable, so that the wavefunction vanishes if two 
particles are at the same point. 

Thus if ^(A^, . . . , A^|Ai, . . . , Ajv; t) denotes the Green's function satisfying the imaginary time 
Schrodinger equation 

then the Green's function G for the Fokker-Planck equation (1.1), with (3 = 2 and W given by (3.1), is 
related to g by 

G(Ar',...,A^|A 1 ,...,A w ;r) = e-^ ^< , ^ ) ))^,...,A^|A 1 ,...,A,;r) (3.5) 

Due to the impenetrability condition of the last paragraph we can suppose that the eigenvalues are 
ordered 

-L/2 < Ai < .. . < X N < L/2, (3.6) 
(and similarly the initial eigenvalues A^), which allows us to write 

p -(w(X 1 ,...,X N )-w(X^\...,X^)) _ TT sin7r(A fc -A 3 )/£ 

For N even (3.7) is anti-periodic under the mapping Xk i— > + L. We then require that g be antipcriodic 
under the same mapping, so by (3.5), G will be periodic. 

It is straightforward to check [6,13] that for N even the Green's function solution of (3.4) satisfying 
the required conditions is 

g(X^ ,...,\ ( £>\\ 1 ,...,\ n ;t) = e^h det[^ 2 ( 7 r(A j - xf)/L- q)U^,..., N (3.8) 

where 

oo 

q:=e -2^r h L^ e 2 (z;q):= £ g (»-l/2) 2 e 2iz(n-l/2) ^ ( g g) 

n— — oo 

(for N odd 62 is to be replaced by 6*3). Thus we have that for N even the required Green's function is 
given by substituting (3.7) and (3.8) in (3.5). 

3.2 The integration technique 

We want to calculate the density-density correlation (2.8) with G given as above and the initial condition 

p(Ai \...,A£>) = J- J] Isin^Af-AfVLI (3.10) 

1N l<j<k<N 

Since we are considering the periodic model the integrations in (2.8) must now be from —L/2 to L/2 and 

p ( i)(A';0) = p (1) (A;r) =N/L = p. 
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The exact value of the normalization C\n is known but will not be needed in our subsequent analysis. 
To calculate (2.8) we will consider the generalized partition function 

N / f L/2 r L/2 \ 

Z(a,b):=l[( / dXfhiX^) d\ l a(\ l ))p(\ { °\...,\ { ° ) )G(\ { °\...,\ { ° ) \\ u ...,\ N] T) 

fJl \J-L/2 J-L/2 ) 

(3.11) 

and use the formula 

S { X^r) = ^ m lo g Z(a, b )_ (3.12) 



i=6=l 



Substituting for p and G we have 

N / r L/2 r L/2 



Z(a,b)=Cmf dX?h(X^)[ dMaiXt)) JJ "S^ " *f) 

1=1 \"'- L / 2 J-L/2 j ]_<j<k<N 

x det[^ 2 ( 7 r(A j - Ai 0) )/i; 9)] j , fe= i,... ; jv det[e 2 ^( fe -( JV+1 )/ 2 )/ I '] j , fe=1 ,..., J v, (3.13) 
where C is independent of A j ^ , A/ , and the Vandermonde determinant formula 

JJ 2sin^(A, - A J )/L = *-^- 1 )/Mct[ e 2 - A ^ fc -( JV + 1 )/ 2 )/ i ] J , fc=1 ,..., Af (3.14) 

l<j<k<N 

has been used. 

To simplify (3.13), we note that since 

det[^ 2 ( 7 r(A j - X^)/L;q)] jik=lt ..., N 

is anti-symmetric in each A; , we can make the replacement 

N 

det[e 2 7 ra J (fe-(W+l)/2)/L ]j fe=i ^ ^ m J-j- e 2,r^0-(iV+l)/2)/i 

in the integrand. The integrations over the A/ can now be performed row-by-row to give 
JL ( r L / 2 

1=1 \"~ L / 2 / l<j<k<N 



Z(a,b) = CN\f[l f dA< >6(A<°>)) [] sgn(Ajj°>-Af )<iet[A j (\M)] j , k=1 ,..., N (3.15) 

1=1 \ J -L/2 j i<j <k < N 

■■= T f L/2 d\a(\)9 2 (n(X xf)/L; q)e ^K)-(N + i)/2)/L ^ (3 lg) 

L J-L/2 

at ions o- 

[1] to give 



where 

rL/2 
-L/2 

The integrations over the a[°^ can now be done using the method of integration over alternate variables 



Hjk : 

Since 



Z(a, b) = CM 2 A C t[H 3k ])'^ N (3.17) 
■ jk ■= f " d\Mb(\W) f L/2 dX^b(X^)A 3 (X^)A k (X^)sgn(X^ A<°>). (3.18) 

J-L/2 J-L/2 



where 

r L/2 j-L/2 



-L/2 J-L/2 

1 S 1 s 

2det[H jk ] SaJX) dct ^ a=6=1 = 2dct[H jk ] 5b(X') dct ^ „ =6=1 ( 3 ' 19 ) 



(3.12) and (3.17) give 



S(A ' A ' ;T) = + 2detk,] fa (A)l(A0 det[ ^ ] L= 6 = 1 - ^ 



3.3 Calculation of the determinants 

From (3.16), and the definition of 62, when a = b = 1 

Aj(X) = e ^(3-(N+i)/2)x/L q u-(N+i)/2f_ (3 21) 



With this value of Aj(X) we have 

r-L/2 



-L/2 

Substituting (3.21) and (3.22) in (3.18) shows that when a = & = 1 



(3.22) 



iJ jfc = hjSjy+i-j-kfi, (3.23a) 

where 



hj ~ _ ^ 2 (j --(iV+l)/2)' (3 236) 



1/ 

7ri(j - (JV + l)/2) 

Hence [i?j/ s ]j J / s =i > ... j Af is only non-zero along the right-to-left diagonal. 

Taking the functional derivative 5/Sb(X') in row j and setting a = b = 1 gives 

I jtk (X') := j ^dX ( - ^(A)A,(A') + A,(A')A fc (A)Jsgn(A' - A) 

= e 2ni(j+k-(N+l))\' /L q (j-(N+l)/2) 2 + (k-(N+l)/2 f ]± ' 1 

7TI 



fc-(JV + l)/2 j-(AT + l)/2, 

(3.24) 

Taking the functional derivative 5/Sa(X) in row j' and setting a = b = 1 gives, for j' ^ j 

Jy, fc (A) := i / L/2 dA(°) / L/2 < ) (0 2 (^(A-A(°))/L; O ) e 2 ^'^ JV + 1 )/ 2 )^A fc (A( O) ) 

^ J-L/2 J-L/2 V 

+0 2 (tt(A - A( 0) )/i;fl)e 2m(fe - (Ar+1)/2)A/i ^(A (0) ))sgn(A(°) - A< 0) ) 
= ^-(^J-e^-^W _ (,' „ *), (3.25) 

while for j' = j we obtain 

: = ie 2 (7r(A - A')/L; g) dA<°> f e 2-0-(Jv+i)/2)A/L Afc(A (o )) _ (j ^ fc) W A ' _ A (o) } 

L J-L/2 V ' 

+ (e**iU-(N+i)/WLA k (\')-(j "k))j j L/2 dX^ sgn(X^ - X')9 2 (7r(X- X^)/L;q) 
= ±9 2 (n(X X r yL . q) ^__L_^ {j „ k) j 

+ ^-(A f +i)/2) 2 e 2 7r4 0A +fe A'-(A + A')(iv + i)/2)/L _ (j ^ fc)U / dA( >sgn(A( ) - A'^MA - A( ))/L;q) 

(3.26) 



From (3.23), with a = b = 1, in each row j, j') the only non-zero term is along the right to left 

diagonal and is given by hj" . Expanding the determinant along this diagonal shows that (3.20) reduces 
to 

5(A',A;r) = -2p 2 + i( £ ^,N + i- S {XVr,N + i-A>) - {j «f, N+1 - j{f) nxcd) ) 

^ j,i' = i J J ' 



+ E ^ ,JV+ y (A ' A/) )- (3-27) 



We note that the above expression is unchanged if we remove the restriction j ^ f in the double sum. 
Doing this, and noting 



N r N 



Ij,N+l-j{X) _ \ - Jj',N+l-j'W 

3 = 1 

(3.27) reads 



1 j,N+i-j{<\ ) _ "J ."t-l-j vv _ 2p 
j=i hj j'=i k 'i 



S(A ,A,r) = 2~(E ^ E ^ J • ( 3 - 28 ) 



The quantities in the above expressions are given explicitly by (3.23b) and (3.24)-(3.26). We see im- 
mediately that the double sum is a Riemann approximation to a definite integral, with the approximation 
becoming exact in the thermodynamic limit. Taking the thermodynamic limit of the first term requires 
an intermediate step: from (3.26) the limiting behaviour of 6 2 is required. 

From the Poisson summation formula 

1 11 00 



^ e -<n-l/2f e 2i(n-l/2)z = ^ (-l) m e - (z+ 

n— — OO TTL— — OO 

Hence with z = 7r(A - A')/L and e = 2tt 2 t/jL 2 



I^^A-AO/i;?) = E (-ire-^'+^) 2 /^ 

^ ' m— — 00 

„ f^) 1/2 e -7(A-V)V2. as i^oo 
V27rr/ 

/OO 
dui e-^Wfh C os 27ra lP (A - A'). 
-OO 



(3.29) 



Using this asymptotic formula, and from the remarks in the paragraph above, we thus have that in 
the thermodyamic limit 

5(A',A;r) = 

/oo /.1/2 
du 1 e- 2T ^ pUl '> 2 /~<cos2iru 1 p(\- A') / du 2 e 2T ^ pUl ^^ cos2iru 2 p(\ — A') 
-co J-l/2 

/co -2t(7tpmi) 2 /7 rl/2 
dui sin27rw 1 p(A - A') / du 2 u 2 e 2r ^ pu ^ ^ sin27ra 2 p(A - A') 
-00 ui J_ 1/2 



P - f 1/2 d«l d U2 ( Ul + U2 )fl e - 2 -(-P) 2 ("?^)/7 + l e -2r(.p) 2 (^- U ?)/ 7 \ 2W ;, ,, V, 

2 J-l/2 7-1/2 W U 2 ) 



e 

/2 

(3.30) 



where to obtain the second term the formula 



/oo />oo 
dA(°)sgn(A(°) - A') / dm e - 2 -(^«i) 2 /7 cos 27rui/>(A<°> - A) 
-00 J —00 

I r°° e -2r(7rpMi) 2 /7 

= / du\ sin27ruip(A — A') 

W-00 U l 



has been used. 
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Expressing the double integrals as single integrals and simplifying gives 



S(X',X;t) :=S{\' -A;r) 

dm e -^pui) 2 /2i cos nUl p(\ _ A') J du 2 e ^(^) 2 /27 cos nU2 p(\ _ A') 

r°° e -T(xpu 1 ) 2 /2"/ r-l 

+ / dui sinTrwip(A-A') / du 2 u 2 e T ^ pu ^ ^ smiru 2 p(\ - A') (3.31) 

Ji u i Jo 

We note that the first term is precisely the density-density correlation for initial conditions with uni- 
tary symmetry given by (2.11), while the second term is precisely that which occurs in the two-point 
distribution function (2.4) (after replacing r by 2r) with (3 = 1 (after replacing r by 2t). The large r, 
A — A' expansion of the above integrals is obtained by expanding the integrands in the neighbourhood of 
ui,u 2 = 1. The behaviour (2.1) with (i = 1 is verified. 

The Fourier transform with respect to A — A' of each of the products in (3.31) have been calculated 
in ref. [7] and [5] respectively. Thus we have 

S(k;r) = §!(k-,T) + S 2 {k-r) (3.32a) 

where 

Si(k;r) := — — exp(-T7r|fc|g max p/7)sinh(r7r|fc|g min p/7) (3.326) 
Tir\k\ 

with 

<7min := min(l, \k\/2np) g max := max(l, \k\/2wp), 

and 

'v 9 \k\ •> r 1 +\ k \/-*p i 

S 2 (k; t) := T-Wl - e~ Tk h) _ U e ^ /2 7 / dfcl e -^,|/#i/7 ; < |fc| < 2^, (3.32c) 



5 2 (fc; r) := — r^e- Tfe ^ sinhTTTplfclAy - ^e Tfe /27 / dfci e -^ P \k\k,h ^ > 2?r/9 _ 

(3.32d) 

Although the small-|fc| behaviour of this expression agrees with (2.13a), we see that there is also a 
singularity at \k\ = 2itp. This corresponds to oscillatory terms in the large-|A — A'| expansion of (3.31) 
and indicates a preference for crystalline ordering with spacing \j2ixp. 
For t, |A — A' | -» (3.31) gives 

S(X' — A;t) ^ / dui e^'/ 2 '' cosWA - A^/r 1 / 2 

Thus, from (2.17), 



Z" 00 2 

u) = / dui e^ 111 ^ 27 costtuiv 
Jo 



/( 

- kr^, ,3.33, 

which is identical to the result (2.20) for initial conditions with unitary symmetry. It is also identical to 
the result found by Simons and Altshulcr [10] when the initial and final distributions are given by (2.2a) 
(i.e. for a perturbation of the eigenvalue p.d.f. with orthogonal symmetry). 



11 



4 Exact solution with symplectic symmetry initial conditions 

The Green's function for the Dyson Brownian motion model for the circular ensemble with (3 = 2 and N 
even is given by substituting (3.7) and (3.8) in (3.5), and is of course valid independent of the particular 
initial conditions. In this section we want to calculate the density-density correlation with the initial 
conditions (2.2b). 

The first step is to take the limits 

Ag^-Ag) j = l,...,N/2 

in the Green's function. We readily find, assuming the ordering (3.6) 

x det[ faWXj X^)/L;q) « ^(A,- - X^)/L;q) } i=1 ,....„ (4.1) 

OA 2fc k = l,...,N/2 

To obtain an explicit expression for the generalized partition function (3.11) we multiply G by 



= f L \ N -rEoh ni< 3 <fc<jv si n^(A fc - Xj)/L 

Ul< j<k <N/2 Sin4 *(>£k ~ X 2j)/ L 



C 4 (N/2) ^k<N/2 



P (xi°\- . . , a®) = 7 ^- z n - 4?)/^ 

This gives 

N ( ,L/2 

Z(a,b):=H[ dA<°>6(A<°>) 

xdet[ ^(Xj-X^/L-q) » ^ 2 (tt(A , - A<°>)/L; q) } j=1 „ det[e 2 ^ k ^ N +^/% 



L/2 N 

dXi a(Xi) 

-L/2 . 



fc = l,...,JV/2 



(4 

where we have used the Vandcrmondc determinant formula (3.14). 

The integrations over Ai , . . . , Ajv can be performed row- by- row while the integrations over A^ , . . . , \$ 
yield a pfaffian structure [1] . We obtain 

(L/2 \ I/ 2 

det[J ^dX^b(X^)Uj k (X^)}j,k=i,...,Nj , (4.3) 

where 

U jk (XM) := A^X^^LjA^) (j - k) 

and Aj is given by (3.16). The density-density correlation can be calculated from (4.3) using (3.20) with 
Hjh replaced by Ujk and — 2p 2 replaced by — p 2 . In fact the calculation is very similar to that presented 
in Section 3.3. 

Thus with a = b = 1 

Ujk = Uj5N+i-j-k,o (4.4a) 

where 

Uj := 2ni(k - j) q ^-(N+i)/2f _ (4 45) 
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so that [Ujk]j,k=i,...,N is non-zero only along the right-to-left diagonal (c.f. (3.23)). Taking the functional 
derivative 5/5a(X) in row j and setting a = b = 1 gives 

Vjk{X ) := ^L e ^H3+k-(N + l))/L ^ _ {N + 1)/2))q 2(k-(N+i)/2 f _ {j „ ^ (4 5) 

while taking the functional derivative S/6b(X') in row j' gives 

W; k (\') = — (k - j') q (k-(N+l)/2) 2 + U'-(N+l)/2) 2 e 2 m \'U'+k-(N+l))/L ( 4 _g) 

for j ^f, and 
Xjk(X, A') 

= ±e^-( N +WV L 2 (n(\ - X')/L; <l)-^M X ') ~ V " k ) 

+ l e 2niX(k-(N+l)/2)/L e 2niX'U-(N+l)/2)/L q U-(N+l)/2f 9 _ y )/£ . g) _ (j fc) (4J) 

for j = j'. 

Since in each row j" the only non-zero element in the denominator of (3.20) (with Hjk 
replaced by Ujk) is along the right-left diagonal and is given by Uj>, we conclude (c.f (3.28)) 

S (\\\-,T) = 1 -(jr Xj ' N+1 - j{X ' X ' ) - ]T Vj.N+i-rWw^N+i-jjX') u „ f) \ (48) 

\ j — i J j j , =\ J J / 

where we have used the fact that 

Vj,N+i-j(\) =2^T W 3'^+i-j'(X') _ 2 p 

i u 3 i i u i 

3 = 1 J 3 =1 3 

The thermodynamic limit of (4.8) is computed analogously to that of (3.28). We find 
5(A',A;r) := S(X' — A; r) 

= ^^(y ^ie- r(7r ^ l)2/27 co S7 rwip(A- A')y du 2 e T(7r " U2)2/2 '' cos7ru 2 p(A - A') 

/■oo /-l T(7r P M 2 ) 2 /27 

+ / duiMie- T(7rpMl) /27 sin7ru 1( o(A- A') / du 2 smiru 2 p(\ - A') (4.9) 

A ' Jo u 2 

This expression is very similar to that for the orthogonal symmetry initial condition given by (3.30). In 

particular, apart from the factor of 1/2, the first term is the density-density correlation (2.12), while the 

second term occurs in the expression (2.4) with f3o — 4 for the two-point distribution. Also, the method 

of obtaining the large-|A — A'|, r expansion is the same as in (3.30). Behaviour corresponding to (2.13b) 

is verified (here it is essential that the large-r limit also be included, for otherwise the leading- 1 A — A'| 

expansion is oscillatory and decays 0(1/|A — A'|), analogous to the t = result [1]). 

The Fourier transform with respect to A — A' of (4.9) is readily calculated. We find 

S(k; t) = \ (&(fc; r) + S 3 (k; r)) (4.10a) 

where Si(k;r) is given by (3.32b) and 

S 3 (k; t) := -J—e-vrpWh sinhTrrplfclAy + lM e -^ 2 /2 7 f J_ & -^ P \k\k^h ^ < |fc| < 2irp, 

7TT|fc| 2tt Ji-\k\/* P fc i 

(4.106) 
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S 3 (k; t) := -1- sinh7TTp|fc|/7 + M e -- fe2 /27 f dkl J_ e -^|fc|fc 1 /7 ) | fc | > 27rp . (4.10c) 

Thus S(k;r) has a logarithmic singularity at |fc| = 27T/5. The first singularity corresponds to the leading 
oscillatory term of the largc-|A — A'| expansion of (4.9) which has period 1/np (the initial spacing between 
the doubly degenerate pairs of eigenvalues), while the second singularity corresponds to the preferred 
crystalline ordering of spacing l/2np between eigenvalues in the final state. 
For r, |A - A'| -> 0, (4.9) gives 

S(X'-X;t) ~ — Vt^ / duie-"'/ 27 cos7rui(A- A')/r 1/2 

1 JO 

+ nKx_y) r°° Uie -„?/2 7 sin7rai(A _ y )/r i/2. 

2r Jo 

Thus, from (2.17) 

1 f°° 2 ^ /*°° 2 

f(v) = - I dui e~ Ul / 2j cosnuiv + - / rfui wie _ " l/ ' 27 sin7ruii; 

2 Jo 2 J 

= «•»> 

which is the arithmetic means of the distributions (2.18) and (2.20). 

5 Equispaced initial conditions 

In this section the equal parameter distribution function (2.1) will be calculated for the solution of the 
Dyson Brownian motion model for the circular ensemble with (3 = 2 and the initial conditions 

xf } = -L/2 + (j - l)/p + u, 0<u< l/p, (5.1) 

which says the eigenvalues are all equispaced (harmonic oscillator type spectrum) at r = 0. 

By integrating the Green's function (3.5) over the initial configuration (5.1) (which as a distribution 
consists of equispaced delta function peaks), we obtain for the corresponding eigenvalue p.d.f. 



p(Ai, . . . , X N ; v, t) = 7V-w/2 e £oT/7 det[-U (tt(A, - (k - l)/p - v)/L; q)] jlk =i,.. 



N 



x Yl amn(Xk- Xj)/L (5.2) 

l<j<k<N 

5.1 Analogy with a classical log-potential system 

Some insight into the expected properties of the two-point equal parameter distribution for (5.2) can be 
obtained from a known identity [14] expressing the determinant in (5.2), averaged over v, in terms of a 
product: 

r-l/p 



pj dv det[^i( 7 r(A j /L - (k - 1)/N + v/L); q 1/N )U=i 

v ' Kj<k<N 
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where 

/ oo \ -(N-l)(N-2)/2 

f N (q) = NN^N-IXN-W (fj (1 _ g 2 fe) \ 

Thus p s is proportional to the Jastrow-type product 

^<3<k<N 

which can be interpreted as the Boltzmann factor of a classical gas with potential energy 

l<j<k<N 

and /? = 1. 

It is well known (see e.g. [14]) that the 6^-function is the solution of the two-dimensional Poisson 
equation in doubly periodic boundary conditions, with period L in the ^-direction and period W in the 
y-direction, and q N — er^ w l L . Furthermore, in the limit L — > oo with W fixed, 

-log 6 1 (Tr\/L;q N ) ~ Vi(A) where Vi(A) := — 7r|A|/W 

while 

— log sin7rA/L ~ Vb(A) where V2(A) := — log |A| 

According to Ornstcin-Zcrnickc direct correlation function theory, the small-fc behaviour of the dimen- 
sionless structure factor 

s(fc) := 1 + - / dx pj 2) {x)e ikx (5.3) 

P J — oo 

of a classical gas is determined by the small-Zc behaviour of the reciprocal of the Fourier transform of the 
potential: 

s(k) ~ — = as k — > 0. 

/?y(fc) 

Now Vi(jfe) ~ 7r/|fc| while V 2 (fc) ~ 0(fc 2 ). Hence 

s(fc) - \k\/n (5.4) 

or equivalently, for A — > oo and r-fixed 

5.2 Exact expression for the equal parameter distribution functions 

To calculate the equal parameter distributions (2.1), we multiply (5.2) by unity in the form of 

r N(N-l)/2 N N/2 det[e 2^k/ N] , =i n 

k = -N/2,...,N/2-l 

to obtain [14], 



r JV(JV-l)/2 



Y[ 2simr(\ k -\ j )/L 

l<j<k<N 



deb[^e^ k+1/2) ^- v)/L 9 3 (irp(^i - u) + in 2 TN(k + l/2)/L 2 T , q" 2 )] j=1 « (5.5) 

Jj k = -N/2,...,N/2-l 
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The advantage of this expression is that with the product expressed as a Vandcrmonde determinant using 
(3.14), it consists of a product of determinants in which the elements have an orthogonality property: 

f L/2 dXj e -2-(fc+l/2)(A 3 -,)/L e 2 7ri (fe' + l/2)(A,-,)/L fe(7rp(A . _ y) + { ^ tN ^ + 1/2 ) / L 2 T , q N " ) = L5 k . k > 
J -L/2 

(5.6) 

This feature allows the distributions (2.1) to be expressed as the determinant 

P(„)(Ai, . . . , X N ; v, t) = &ct[K (Xj, X k ; v, T)]j jk =i,...,N (5.7a) 

where 

JV/2-1 

K(Xj,X k ;u,T) := ]T -e 2 ^ +1 ^-^' L ^ 3 {-Kp{X k - v) + in 2 TN(p + l/2)/L 2 7 ; e"^) (5.76) 

p=-N/2 

The thermodynamic limit follows by noting that this expression is a Riemann sum approximation to a 
definite integral, which becomes exact in the thermodynamic limit, thus giving 

.1/2 
1/2 

From (5.7a) and (5.8), the truncated two-point distribution averaged over v is given by 

,1/2 ,1/2 

ds 

1/2 J-l/2 

xd 3 (irp(X - A') + iir 2 T P 2 (s - s')/ T , e^ 2 " 2 ^) (5.9) 



,1/2 

K(Xj,X k ;u,T)=p / d S e- 2m P^- x ^ 8 9 3 (irp(X k - v) +nr 2 T P 2 S / T ,e- 27T2p2T ^) (5.8) 

J-l/2 

a) and (5.8), the truncated two-point distribution average 
P J 2) (X-X';t) - -p 2 [ 1/2 ds [ 1/2 ds ' e -^-y)p(s-s') 

J-l/2 J-l/2 



From this expression, it is straightforward to compute the small-|fc| behaviour of the structure factor (5.3) 
(only the n = 0, ±1 terms in the series definition of 6*3 contribute). We find 

*(*) ~ M(l +e - 2 - 2 -Pl fe l/7) ( 5 . 10 ) 

For fixed-r the term in brackets is to leading order equal to 2, so the result (5.4) is verified. 
5.3 Asymptotics of the spacing distribution 

The distribution p(s) of the nearest neighbour spacing averaged over v can be computed by differentiation 
of the probability h(s) that an interval of length s averaged over v is free from eigenvalues: 

f (-) = -|jM«)- 

The leading asymptotics of h(s) can be predicted by adapting an electrostatic argument of Dyson [15]. 

In Section 5.1 we have noted that the eigenvalue p.d.f (5.2) averaged over v can be interpreted as the 
Boltzmann factor of a classical gas with potential which to leading order in large-A behaves as 

Vi(A) := -7r|A|/W, W = 2ttt P /j. (5.11) 

Now (5.11) is proportional to the potential of the Coulomb repulsion between sheets of charge (i.e. the 
one-dimensional Coulomb potential), and the corresponding classical gas is the one-component one- 
dimensional plasma. The basis of the electrostatic argument is that in the macroscopic hole size limit 
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s — > oo, the plasma behaves like a perfect conductor. Thus the total charge — ps inside the hole due 
to the uniform background will be exactly cancelled by two delta function charge distributions at each 
boundary of the hole of charge ps/2. 
According to the hypothesis of Dyson 

h(a) ~ e -(£M«)+tM«)+tfsM) (5.11) 
where U\(s), ^(s) and 1/3(3) are given by 

C/i(s):=V^f J o dx J Q dy\x-y\ (5.12a) 

U2(s):=p 2 ^(J S dx\x\+ J S dy\ S -y\) (5.126) 

««•> <v*> 

and represent the electrostatic energy of the background-background, background-boundary and boundary- 
boundary interactions respectively. Evaluating (5.12a) and (5.12b) and substituting in (5.11) gives 

h( S ) ~ e -(ps) 3 7/24p\ ( 5J3) 

This is to be contrasted to the asymptotics of h(s) with r — > 00 (see e.g. [1]): 

h( S ) ~ e -(^) 2 /8. 

6 Multiple parameters 

Beenakker and Rcjaci [7] considered the multiple-parameter Fokker-Planck equation 

with £ given by (1.1) and W by (1.2), as a model for the statistical description of the dispersion relation 
of a d-dimensional crystalline lattice. With the initial condition (2.7) extended to f — 0, it is easy to 
see that the Green's function solution of (6.1) with = 2, N odd and W given by (3.1) is obtained by 
making the replacement 

ie 3 (7r(A, - \ k )/L- q) » - d J2 WA,- - A fc )/L; g„), ^ := e" 2 ^/^ (6.2) 

Furthermore, explicit calculation shows that the density-density correlations and equal-time corre- 
lations for the multiple-parameter equation are the same as for the single-parameter equation (with 
corresponding initial eigenvalue p.d.f.'s) provided the replacements 

e±T Vu 2 T/2 7 „ ^_ e -»V« 2 W27^ ( 6 .3) 

are made in each of the integrands of (2.11), (3.32) and (4.9), while in (5.7) the replacement 

1 d 

3 (7rp(X k -v) + iw 2 rsp 2 / r , e-^e^h) 3 (np(X k - v) + iw%sp 2 / 7fl ; e - 2 -VW7.) (6.4) 
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is required. 

As an example, the multiple parameter generalization of the density-density correlation (2.11) for 
initial conditions with unitary symmetry is 

S{\-\';t)= p 2 J™ d Ul ^I^ e --V«?W27^ cos7ruip(A-A') 

* Jo ^ ^ e " 7r2p2 "' T '' /27 ^ cos7rw 2 p(A- A') (6.5) 

For large- 1 A — A'| and t m (fi = 1, . . . , d) (with each r M of the same order), the asymptotic expansion of 
(6.5) is obtained by writing 

£ =exp(-l£ ^(^) 2 ) 



By the assumption on the relative order of each r^, we see the second term on the r.h.s. is O(l) and 
can so be ignored for purposes of computing the leading order asymptotic behaviour. Thus we again find 
the formula (2.10) with = 2, provided the replacement 



Jl^IyjjL (6 . 6) 

2 7 d^2 la { ' 



2 7 d £j 27„ 

is made. This same prescription gives the asymptotic expansion of the density-density correlation with 
multiple parameters for orthogonal and symplectic symmetry from the corresponding expansion for the 
single parameter case obtained in the text. These finding are in agreement with those in [7] (however the 
requirement that the expansions are only valid for each t m of the same order was not made explicit in 
[7])- 

Similarly, we find the distributions (2.17), (2.19), (3.33) and (4.11) are unchanged, provided the 
replacement 

(S ,/2 ^"i£© 1/v ^ <«> 

is made. 



Note added: Since the completion of this work, I have been informed by J. A. Zuk that the results of ref. 
[3] include the calculation of the density-density correlations given in Sections 3 and 4. However these 
results, which are derived using the supersymmetry method, are presented as complicated triple integrals, 
whereas our results (3.30) and (4.9) are products of simple single integrals. We have not been able to 
show directly that the different forms are equivalent. Instead Zuk [16] is repeating the calculation of [3], 
with a different parametrization [17] which is known to reproduce the equal parameter results (2.4)-(2.6), 
in a bid to derive (3.30) and (4.9) using the supersymmetry method. He also points out that taking the 
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Laplace transform of (3.32) in the variable t := (7rp) 2 r/27 gives a result entirely in terms of elementary 
functions: 

/OO r-OG 
dxe ikx dte- st p- 2 S(X;T), (x := irpX) 
-oo JO 

( s log fc 2 log(l + W). 1*1 < 2 

2\k\(s-k>)\ S logJ|±i^-^log|±l, |fc|>2 
We also note that the paper [18] has appeared which states a result (eq. (47)) equivalent to our (5.9). 
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